arXiv: 1503.03995v2 [cond-mat.mes-hall] 2 Sep 2016 


Simulation of supersymmetric quantum mechanics in a Cooper-pair box shunted by a 

Josephson rhombus 

Jascha Ulrich, Daniel Otten, and Fabian Hassler 
JARA-Institute for Quantum Information, RWTH Aachen, D-52074, Germany 

(Dated: November 2015) 

Supersymmetries in quantum mechanics offer a way to obtain degeneracies in the excitation 
spectrum which do not originate from selection rules. The mechanism behind the degeneracies is 
the same as the one that leads to the miraculous cancellations of divergences in supersymmetric field 
theories found in the high energy physics context. There is up to now no realistic proposal of non- 
integrable systems that show level degeneracies due to a supersymmetric structure. Here, we propose 
an implementation of a quantum-mechanical supersymmetry in a Cooper-pair box shunted by a 
Josephson junction rhombus which is effectively 7r-periodic in the superconducting phase difference. 

For a characteristic ratio between the strength of the 2n- and the 7r-periodic junction, we find a 
two-fold degeneracy of all the energy levels all the way from the weak junction/charge qubit limit to 
the strong junction/transmon regime. We provide explicit experimental values for the parameters 
of the system and show that tuning in and out of the supersymmetric point is easily achieved by 
varying an external gate voltage. We furthermore discuss a microwave experiment to detect the 
supersymmetry and conclude that it can indeed be implemented with currently existing Josephson 
junction technology. 
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The macroscopic quantum mechanics of supercon¬ 
ducting circuits has allowed the experimental simula¬ 
tion of many complex quantum phenomena such as 
phase transitions , 1 quantum spins , 2 or dynamics in open 
systems . 3 Theoretically, the quantum simulation of in¬ 
tricate subjects such as Hawking radiation 4 and lattice 
gauge theories 5-7 has been proposed. In the plethora of 
phenomena that can be simulated with the help of su¬ 
perconducting circuits , 8,9 degeneracies due to quantum- 
mechanical supersymmetries have notably been absent. 
Typically, degeneracies in the spectrum arise when the 
Hamiltonian commutes with all group elements of a non- 
Abelian symmetry which translates into selection rules 
demanding vanishing off-diagonal and equal diagonal ma¬ 
trix elements of the Hamiltonian within the same irre¬ 
ducible representation . 10 The degeneracy of the states 
thus always follows from the dimension of the representa¬ 
tion. Supersymmetry on the other hand does not simply 
forbid different states to couple but it makes sure that in 
each order of perturbation theory there is always a pair 
of terms canceling each other . 11 ’ 12 

It is intriguing that the degeneracies of supersymmetric 
quantum mechanics occur by the same mechanism 11 that 
leads to a miraculous cancellation of divergences in su¬ 
persymmetric field theories and makes supersymmetries 
an important tool of particle physics . 13 In trivial cases 
like the free particle 14 or the Jaynes-Cummings model , 15 
however, the supersymmetric structure is irrelevant since 
the spectrum is exactly solvable. In order to deepen the 
connection to the ideas in the high-energy context, it is 
thus of vital importance to propose a non-integrable sys¬ 
tem where the level degeneracy can be solely understood 
by its supersymmetric structure. 

In this paper, we show that shunting a Cooper-pair 
box with a Josephson junction rhombus simulates a 
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FIG. 1. Setup simulating the Hamiltonian Eq. (1). A Joseph¬ 
son junction (crossed box) with Josephson energy Ej\ and 
capacitance Ci couples the Cooper-pair box, a superconduct¬ 
ing island with superconducting phase <j>, to a ground su¬ 
perconductor (diagonally striped) with phase <fo = 0. A 
Josephson rhombus (depicted as a rhombus) provides an addi¬ 
tional shunt to the ground which generates a cos(2 rf>) Joseph¬ 
son coupling of strength Ej 2 when the rhombus is threaded 
by a flux = <E>o/2. We assume that an additional flux 
<&jr = —$_r/2 threads the loop between the standard junc¬ 
tion and the rhombus, which requires a second flux line that 
can be controlled separately. A capacitance C g couples the 
system to a transmission line at the voltage V(t). 


highly nontrivial supersymmetry that can be experimen¬ 
tally realized with today’s Josephson junction technology. 
This proposal combines ideas of implementing supersym¬ 
metry in purely bosonic systems 16 with the quantum- 
mechanical supersymmetry that has recently been pro¬ 
posed for superconductors hosting fermionic Majorana 
bound states . 12 The only nonstandard component of 
our proposal is the Josephson rhombus. The Joseph¬ 
son rhombus is a junction between two superconduc¬ 
tors that allows only tunneling of pairs of Cooper pairs. 
Consequently, its current-phase relation is Tr-periodic. 1 ' 
















2 


Josehson rhombi have previously been proposed as build¬ 
ing blocks for topologically protected qubits 18,19 that 
have been shown to be experimentally feasible. 20 Ad¬ 
ditionally, they have been employed for the experimen¬ 
tal realization of qubits based on the Cooper-pair par¬ 
ity. 21 Up to now, theoretical studies on the Josephson 
rhombi have mainly been focused on the semi-classical 
regime. 22,23 For our setup, we study a Josephson rhom¬ 
bus in the charging limit where it is adiabatically coupled 
to the superconducting island and generates a 7r-periodic 
Josephson coupling of a specific cos(2 <j>) form. A super- 
symmetry is then obtained for a characteristic ratio be¬ 
tween the strength of the conventional and the 7r-periodic 
Josephson junction, see below. 

Our system of interest is depicted in Fig. 1. It is an ex¬ 
tension of the conventional Cooper-pair box, 24,25 which 
consists of a superconducting island with superconduct¬ 
ing phase (f> and a ground superconductor with phase 
</>o = 0 which are coupled by a Josephson junction with 
Josephson energy Eji and capacitance C±. A capaci¬ 
tance C g couples the system to a transmission line biased 
at a DC voltage V(t) = V g . We add an additional shunt 
to the ground through a Josephson rhombus with capac¬ 
itance Co, which, as we will discuss in more detail below, 
generates a 7r-periodic Josephson energy — Ej 2 cos(2</>) 
when threaded by a flux <!># = <h 0 /2, where $o = h/ 2e 
is the superconducting flux quantum and Ej 2 the effec¬ 
tive junction energy. Taking into account an additional 
flux Qjr = — &r/2 in the loop between the conventional 
junction and the rhombus, we obtain the effective low- 
energy Hamiltonian 

H eS = 4E Cs (n - n g ) 2 - E J1 cos (f> - Ej 2 cos(2</>), (1) 

where n = —id/d(f) is the number of Cooper pairs on the 
island, Ec s = e 2 /2Cs with Cs = Ci+C g +C 0 is the total 
charging energy of the island, and n g = C g V g /2e is the 
induced offset charge in units of 2e. The Hamiltonian (1) 
does not admit an analytic solution. Its only symmetry is 
the operation K : i—>- —0 at the point n g = 0 which due 
to its Abelian nature does not lead to any degeneracy. 
However, as we will show below, for a specific ratio of 
the energy scales all excited levels are degenerate due to 
a supersymmetry. 

In the simplest setting, a Hamiltonian Hq is called 
supersymmetric when there exists a Hermitian involu¬ 
tion K with K 2 = 1 that commutes with Hq and a 
Hermitian supercharge Q which anticommutes with K 
and factorizes the Hamiltonian Hq = Q 2 . 26,2 ' The sec¬ 
tors of Hq are then characterized by K according to 
Hq = P+H q P + + P-HqP_ with P± = (1 ± K)/2 and 
are intertwined through the relation P±Q = QP^i which 
guarantees that to each eigenstate |a) to energy E a > 0 
in one sector there is a partner state (Q/y/H,)\a) to the 
same energy in the other sector. To see how this relates 
to our system, let us introduce the supercharge Q and 
the involution K according to 12 

Q = 2\jEcv (n + msin</>) (— l) n , K</>K = —</>, (2) 


(a) (6) 



FIG. 2. (a) Circuit of the Josephson rhombus consisting of a 
loop interrupted by four Josephson junctions with Joseph¬ 
son energies Ej 0 and capacitances C 0 . When the loop is 
threaded by a flux of = d>o/2, the transport of single 
Cooper-pairs through the device is suppressed by destructive 
interference between the tunneling events through the left and 
the right arm of the rhombus, rendering the rhombus eigenen- 
ergies 7r-periodic in the fixed phase difference (j> between the 
top island and the ground superconductor. ( b ) Choice of the 
gauge-invariant phase differences across the links with the ar¬ 
rows indicating their orientation. The phases 9j are conjugate 
to the charges rij = —id/ddj of the superconducting islands 
(gray squares) in the two arms of the rhombus. The phases 
add up to the reduced flux <f>R = 27t$_b/$o- 

where a is a free parameter. Note that the supercharge Q 
is Hermitian since the addition/removal of a Cooper- 
pair anticommutes with the Cooper-pair parity (—1) of 
the island. We then find that with a = \JEj2/2Ec-s and 
up to irrelevant constants, the supercharge Q squares to 
the Hamiltonian Eq. (1) at the point n g = 0 and 

Eji = y/8 E J2 Ec s , (3) 

where the system Eq. (1) is supersymmetric. The exotic 
feature of the supersymmetry to note here is the preser¬ 
vation of the degeneracy of the energy levels all the way 
from the charge qubit regime with a -C 1 to the transmon 
regime a > 1 as long as (3) is fulfilled. The “hidden” 
character of this degeneracy is underlined by the highly 
nonlocal form of the supercharge (2). 

From the above, we see that an (effective) Joseph¬ 
son junction with a 7r-periodic Josephson energy of the 
form — Ej 2 cos(2</>) is crucial for supersymmetry. Such 
a circuit element is provided by the Josephson rhombus 
shown in Fig. 2(a). 1 ' It is a two arm Cooper-pair inter¬ 
ferometer connecting the superconducting island to the 
ground in which single Cooper-pairs tunneling through 
the left and right arm of the rhombus interfere destruc¬ 
tively due to a magnetic flux = < Fo/2. Each arm con¬ 
tains two Josephson junctions connected in series with 
Josephson energy Eco and capacitance C 0 . We show 
that in the charging regime ry = Ej <> /2Ec<> < 1, the 
ground state energy of the rhombus is well approximated 
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FIG. 3. Energy spectra of the full Hamiltonian Eq. (4) as a function of the dimensionless offset charge n g = C g V g /2e for 
a = 2,r/ = 0.7 in (a) and a = l,r/ = 1 in (c). For a given r/ = Ej g /2Ec<>, we numerically determine the strength Ej 2 
of the 7r-periodic component of the ground state energy eo (</>) of the rhombus. The value of the charging energy is then 
obtained from the relation Ec s = Ej 2 /2a 2 . The strength of the Josephson coupling Ej\ is fixed at a value AEc s a which 
corresponds to the supersymmetric point, cf. (3). It can be seen in (a) that all the excited levels cross at n g = 0, confirming 
the validity of the effective supersymmetric model (1) beyond the perturbative regime rj -C 1. As shown in (c) the former 
supersymmetric level crossings at n g = 0 turn into avoided crossings at a slightly elevated r). This signals the breakdown of the 
rhonrbus/island decoupling and thus restricts the mapping of the full Hamiltonian to the supersymmetric model. Additionally, 
the first excited rhombus level can be seen as a horizontal line around E/Eco = 0. ( b ) Plot of the power loss P sc at frequency uj 
of the transmission line coupled to the system with parameters given in (a) which corresponds to the experimental parameters 
(Ci +C g )/Co = 80, Ej i = 0.1 Eco, and Eco = 500 /reV, cf. (6). The power loss is measured in units of the injected intensity Po 
and a dimensionless coupling constant ac- Here, we assume that the system once exited relaxes fast with a rate F = 0.5/reV/fi 
into degrees of freedom different from the transmission line such that the resonance condition is indicated by a dip of size P sc in 
the transmitted intensity. Since the ground state is nearly insensitive to changes in n g , the power loss can be directly compared 
to the spectrum shown in (a). 


by eo(<W ~ — Ej 2 cos( 2 (/>), where </> is the (fixed) phase dif¬ 
ference between island and ground. Furthermore, we ar¬ 
gue that the weak coupling to the island permits an adia¬ 
batic decoupling leading to the effective Hamiltonian (1). 
To this end, let us denote the Cooper-pair number of 
the superconducting islands in the left and right arm by 
rii, n 2 and choose the gauge-invariant phase differences 
across the junctions as indicated in Fig. 2(6). Taking into 
account the additional flux <&j R = —$ 0 /4 from Fig. 1, 
the Hamiltonian of the full system assumes the form 

H = 4Ec s {n — n g ) 2 — Ej\ cos</> + H 0 , (4) 

which corresponds to the effective model Eq. (1) with 
—Ej 2 cos(2 </>) replaced by the rhombus Hamiltonian 

H 0 = 2E C o{nl+nl) + V (> (e 1 ,e 2 ,^. ( 5 ) 

Here, Eq o = e 2 /2 C 0 is a charging energy, the phases 6±, 
0 2 are conjugate to m = —id/ddi, n 2 = — id/dd 2 , and 
the potential V 0 reads 

Ho = -Ej* Y, [cos(| + 0j) + cos(| + (-1)^ - %)], 
1 = 1,2 

where cj) R = 2 tt<& r /<& 0 is the reduced flux. For fixed 
</>, the Hamiltonian H 0 possesses instantaneous eigen¬ 
states | n;<fi) with eigenvalues e n {4>). The destructive in¬ 
terference of single Cooper-pair tunneling is expressed 
by the fact that exchanging the two tunneling paths and 
advancing (/> by it is a symmetry of the system , 28 de¬ 
manding 7 r-periodicity of e„ (</>). At half a flux quantum, 


time-reversal is an additional symmetry demanding an 
even e ra (</>). Consequently, the ground state energy eo (</>) 
must be of the form eo {<j>) = —J2 n Ej2nCOs(2n^i). By 
perturbation theory in 77 , we find that the desired 7r- 
periodic component Ej 2 /Ec<, = 7 ? 7 4 /4 dominates, with 
Eja/Eco = -68687t/ 736864 and E J2m /E C o oc rf m . In 
the following, we are interested in the regime 7/ < 1 
and thus we determine Ej 2: Ej 4 numerically from eo (</>) 
whenever needed. We find that Ej 2 stays at least an or¬ 
der of magnitude larger than Ej\ up to 77 « 1 . Since the 
above form of the rhombus energies is due to symmetry, 
differing Josephson couplings E J<> L / R = Ej 0 ( 1 ± S/2) 
in the left and right rhombus arm will in general induce 
a finite 27r-periodic Josephson coupling whose strength 
scales perturbatively as 4 t/ 2 i5 Eco■ Comparison with Ej 2 
yields that the effects of asymmetry are negligible for 
7/ 2 » 5 . 

Projecting the Hamiltonian Eq. (4) onto the instan¬ 
taneous rhombus ground state | 0 ; <^) and using eo (<j>) ~ 
—Ej 2 cos(2 (/>) leads by standard methods 29 to the ef¬ 
fective Hamiltonian iJ ac j = H e g + ^Ec-^ X7>o | A-nO | , 
where A nm = i(n; <f>) is the induced vector po¬ 

tential describing the nonadiabatic corrections. In de¬ 
riving H a g, we have used that the term Aqq vanishes 
since the states |77i; <fi) can be chosen real. Due to time- 
reversal and rhombus symmetry, |A„o| is even in cf) and 
thus does not couple the supersymmetric partners at 
n g = 0. Since the gap to the next pair of supersymmet¬ 
ric states is at least of order Ec s , the effects of |A n0 | 2 = 
|(n; 01 <9<t>TT 0 10 ; <)>}| 2 /[e ra (<)>) —e 0 (<)>)] 2 ex t / 2 are negligible for 
















4 


ijC 1. The coupling to the excited rhombus levels that 
we do not take into account in the projected Hamiltonian 
i? a d is suppressed by even higher orders in 77 . 

The supersymmetry becomes trivial for a —> 0 where 
one recovers the supersymmetry of the free particle . 14 We 
therefore aim for the most interesting regime of a ss 1 , 
where all the terms in the Hamiltonian (1) are of the same 
order. While the adiabatic decoupling of the rhombus is 
most robust for large scale separation 77 -C 1 between the 
rhombus and the island, our perturbative results for Ej 2 
show that this also implies a oc rj 2 -C 1. The regime 
a ~ 1 thus requires moderately large 77 for which we nu¬ 
merically show that the adiabatic decoupling still works. 
Figure 3(a) shows the numerical results for the spectrum 
of the full Hamiltonian as a function of the offset charge 
n g for a = 2 and 77 = 0.7. We highlight that all excited 
levels, including the levels higher in energy not visible in 
Fig. 3(a), become doubly degenerate as n g approaches 
zero. This degeneracy of all excited states in complete 
absence of selection rules gives a clear signature of super- 
symmetry. Figure 3(c) shows the high sensitivity of the 
supersymmetry to the choice of 77 . For a = 1 and 77 = 1, 
supersymmetry at n g = 0 is clearly destroyed by non- 
adiabatic corrections in the Hamiltonian H a d. Moreover, 
the first excited level of the rhombus showing up as a 
horizontal line in the upper region of Fig. 3(c) is visible. 

The spectrum of the system can be read out with the 
help of the transmission line coupled to the island by in¬ 
jecting a voltage V(t) = V g + K,cos (wt) with the AC 
amplitude V u at the frequency u. For a transmission line 
characterized by the admittance Yti, this corresponds to 
an average injected intensity P 0 = \Y t \V%. According to 
the Hamiltonian (4), the AC voltage drives transitions 
in the system through the coupling He = gcY(t)n with 
gc = 2eC g /C'z- Due to the coupling to the charge n, the 
driving is most effective for a > 1. We assume that the 
system relaxes dominantly into channels different from 
the transmission line with a rate T r a b s , where r a b s 
is the rate of photon absorption. Measuring in trans¬ 
mission, the absorption of photons is then signaled by a 
reduced transmitted intensity Pt with respect to the in¬ 
coming intensity P 0 . With T r a b s , photons of energy 
hio exclusively drive transitions from the ground state to 
excited states and the scattered intensity P sc = P 0 — P t 
follows as P sc = HujT a b s . For the calculation of r a b s , we 
employ Fermi’s golden rule. We find 


P sc = 4acPo ^2 


n> 0 


/i 2 wr|(7r|7<90|O)| 2 
(E n -Eo- M 2 + vr 2 ’ 


( 6 ) 


denoting the eigenstates to energy E n of the full Hamil¬ 
tonian (4) by |ti}; here, the dimensionless constant ac 
is given by ac = gc/KY t \ = 2nC g /C^Z 0 Y t \ with the su¬ 
perconducting impedance quantum Zq = /i/4e 2 ss 1 kf2. 
The condition r a b s -C T translates into gcVui "C tuj- The 
validity of Fermi’s golden rule for the rate calculation on 
the other hand demands acoj -C T. 

The system with a = 2 presented in Fig. 3(a) corre¬ 
sponds to the experimental parameters 77 = 0.7, (Ci + 


C g )/C 0 = 80 and Ej\ = O.lE'co- In Fig. 3(6), we show 
the numerical results for the scattered intensity P sc as 
a function of the offset charge n g and the radiation fre¬ 
quency ui with the system parameters stated above. We 
assume the experimental scale Eco = 500 geV. As visi¬ 
ble from Fig. 3(a), the ground state is almost insensitive 
to changes of n g at a = 2 such that the scattered inten¬ 
sity corresponds directly to the energy spectrum of the 
system. The scattered intensity is strongest for the first 
degenerate pair of levels which reflects the fact that states 
higher in energy show a behavior closer to the charging 
regime than the low-energy states. Note that the whole 
spectrum and the level crossings of all excited states at 
the supersymmetric point (3) with n g = 0 can be clearly 
observed. 

As a last point, let us comment on the susceptibil¬ 
ity to imperfections in design. The above analysis was 
based on the rhombus symmetry which is violated both 
by stray offset charges or parameter variations within the 
offset arms as parametrized by 5. As we explain in more 
detail in the Appendix, the susceptibility to stray offset 
charges can in practice be reduced by adding inductive 
shunts to the ground within the rhombus arms . 19,21 For 
the system parameters chosen above, numerical checks 
show that symmetry violations S are tolerable up to a 
few percent. The most demanding experimental require¬ 
ments are thus the reproducible parameter of the rhom¬ 
bus junctions. In contrast, deviations in Ci/C<> can be 
accounted for by the Josephson coupling Ej 1 which can 
be tuned easily. 

To conclude, we have shown that a quantum- 
mechanical supersymmetry arises in a Cooper-pair box 
when it is shunted by a Josephson rhombus operated in 
the charging regime, where it is weakly coupled to the 
Cooper-pair box and generates an effective 7 r-periodic 
Josephson coupling of the form cos(2 cf>). The supersym¬ 
metry is nontrivial since there are no selection rules and 
an analytic solution is impossible, but yet, it guaran¬ 
tees an exact degeneracy of all excited levels. We have 
shown that the supersymmetry can be detected through 
microwave spectroscopy and tuning in and out of the su¬ 
persymmetric point is easily possible by tuning the gate 
voltage. We have proposed realistic device parameters, 
paving the way for an experimental exploration of the 
exotic degeneracies brought by supersymmetries in the 
near future. 
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well as useful discussions with Christoph Ohm. 


Appendix A: Device with inductive shunts 


The simple rhombus design that we have presented in 
the main text for clarity of the discussion suffers from 
the fact that the charge configuration on the capaci¬ 
tors of the central islands in the rhombus is strongly 
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inductive shunts remain good variables. This means that 
El should be a weak perturbation with El ■C Ej 0 , Ec<>- 
We can therefore follow the ideas of Refs. 30 and 32 and 
transform to a basis of Bloch waves. Solving the 27 t- 
periodic part of the rhombus perturbatively in 77 and 
projecting on the lowest band of the rhombus yields a 
Hamiltonian of the form 



FIG. 4. Rhombus design as proposed theoretically in Ref. 19 
which includes inductances L that shunt the islands in the 
rhombus arms and reduce the susceptibility to stray offset 
charges coupling capacitively to the islands. Modeling the 
offset charges through fluctuating voltage sources that charge 
capacitors C s coupled to the islands yields the relation ni 11 = 
C„V s U) /2e for the stray offset charge n^\ Apart from the 
change in rhombus design, the proposed setup remains the 
same as in Fig. 1. Importantly, it is sufficient to add inductive 
shunts only to the lower part of the rhombus arms such that 
no inductive coupling to the main island is generated. 


susceptible to fluctuations of stray offset charges that 
couple capacitively to the islands. They can be mod¬ 
eled, e.g., through voltage sources that charge 

capacitors C s coupled to the rhombus islands, yield¬ 
ing = CsVs 1 ^^/2e for the stray offset charge 

; see Fig. 4. These stray offset charges are not con¬ 
trolled and fluctuate in general independently, destroy¬ 
ing the symmetry of the rhombus arms and lifting the 
7 r-periodicity of the rhombus. As was argued theoret¬ 
ically 19 and tested experimentally , 21 this can be reme¬ 
died by adding inductive shunts of strength L to the 
lower parts of the rhombus arms. In the classical case 
and ignoring the Josephson junctions, adding inductive 
shunts to the ground in the central islands of the rhom¬ 
bus reduces this susceptibility by a factor of w 2 /w 2 c , 
where ui is the frequency of offset charge fluctuations and 
u>lg = 1/V^LCo = y/4EcoEL/his the plasma frequency 
of charge oscillations in the resulting LC resonator with 
the inductive energy El = (&o/2n) 2 /L. As was shown in 
Ref. 30, this property carries over to the quantum case. 
As a consequence of the inductive shunts, the rhombus 
does not couple any more directly to but only 

to hs ^ (2) . The associated noise power changes from the 
1/w form typical for fluctuations of ni 1 ^*' 2 ^ 31 to a much 
more benign noise power proportional to ui. 

As we have discussed above, the inclusion of inductive 
shunts reduces the sensitivity to offset charge fluctua¬ 
tions, but it should of course also preserve the behavior 
of the rhombus as an effective cos(2 (f>) element. This 
means that in our setup, we cannot make El arbitrarily 
large since our treatment required that both the band 
indices and the Bloch momenta qi associated with the 
translational symmetry in the variables di in absence of 


+ g(qj) cos(20) + (-l) J %j)sin(0) j, (Al) 

with periodic functions /, g , h with period 1. The ab¬ 
sence of a coupling between qi and q 2 reflects the fact that 
the phases 9j are not coupled by the rhombus Hamilto¬ 
nian. The presence of a coupling to cos(2 (f>) [sin(0)] that 
is even [odd] under the exchange q\ F 4 q 2 reflects the 
symmetry under exchange of the rhombus arms and si¬ 
multaneous advance of 4> by 7 r that we discussed in the 
main text. The additional time-reversal symmetry at half 
a flux quantum forbids a coupling to cos(</>) or sin( 2 </>) 
and requires that the functions /, g , h possess quadratic 
expansions around q = 0 which we find to be of the form 

f(q)/Eco = 2 q 2 {l + 0 (r ] 2 : q 2 )) + 0 {g 4 ), 

g(q)/Eco = 1 + ^-q 2 + OfaV)), (A2) 

h{q)/E C o = -877 V (l + 0 (rj 2 ,q 2 )). 

Obviously, for <71 = <72 — 0 and El = 0, the Hamiltonian 
If' reproduces the ground state energy eo (4>) that was 
given in the main text. For finite El, the former rhom¬ 
bus eigenstates with sharp Bloch momenta q\ , <72 are re¬ 
placed by eigenstates of the Hamiltonian (Al). Assuming 
a symmetric state of the rhombus arms, the ground state 
energy t' 0 (<(>) of can still be expanded in the form 
e'o(<t>) = -jy n E J 2 n cos(2 ncj)). Importantly, E l < E C o 
implies q 2 ~ El /4 Eco -C 1 which guarantees that the 

term E ' J2 remains dominant in the expansion of £q 
N umerical checks show that this property persists be¬ 
yond the perturbative regime in 77 and remains valid also 
for 77 = 0.7. The inductive shunts thus manage to re¬ 
duce the susceptibility to charge noise without spoiling 
the generation of an effective cos( 2 ^>) potential by the 
rhombus. 

Finally, let us comment on the adiabatic decoupling 
of the island and the rhombus for finite El- The level 
spacing of the eigenstates of the Hamiltonian (Al) is of 
the order of the plasma energy Hlolc- In view of the adi¬ 
abatic decoupling, it is desirable to have Hlolc ^ A^SUSY, 
i.e., El > e ‘susy/E co, where Asusy is the energy scale 
of the rhombus. We note that for a = 2 and 77 = 0.7, 
we have Agygy ~ Eji = Aco /10 such that choosing 
A L > Agugy/Aco = Aco/100 and E L -C E C o,E Jo is 
easily possible. Following Ref. 21, such an inductance can 
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in practice be implemented through a chain of Josephson junctions. 
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